Chapter 1 

Theory of deflagration and fronts of 
tunneling in molecular magnets 

D. A. Garanin 



Abstract Decay of metastable states in molecular magnets (MM) leads to 
energy release that results in temperature increase that, in turn, positively 
affects the decay rate. This is the mechanism of recently discovered magnetic 
deflagration that is similar to regular chemical burning and can propagate in 
a form of burning fronts in long MM crystals. Near spin-tunneling resonances 
the decay rate is also affected by the dipolar field (self-consistent with the 
switching magnetization) that can block or unblock tunneling. There are non- 
thermal fronts of tunneling in which the magnetization adjusts in such a way 
that the system is on resonance within the front core, so that the tunneling 
front can propagate. In general, both dipolar field and temperature control 
fronts of quantum deflagration. The front speed can reach sonic values if a 
strong transverse field is applied to boost tunneling. 



1.1 Introduction 



Deflagration or burning is decay of metastable states accelerated by the 
temperature rise due to the energy release in this process [lj [2] . In most 
cases the decay rate has the Arrhenius temperature dependence, r = 
Jb exp [—11/ (ksT)], where U is the energy barrier. Because of the very strong 
positive feedback, burning can have a form of a thermal runaway: almost 
undetectably slow relaxation at the beginning followed by an explosion-like 
relaxation at the end (explosions at ammunition-storage sites, Bhopal disas- 
ter, etc.). In other cases there is a burning front propagating with a constant 
speed away from the ignition point. These fronts are driven by the heat con- 
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duction from the hot burned region to the cold unburned region before the 
front. Burning of a sheet of paper is a good example of a deflagration front. 

Molecular magnets (MM), of which the most famous is Mn^Ac (3], are 
burnable materials because of their bistability resulting from a strong uni- 
axial anisotropy that creates an energy barrier One can make magnetic 
state metastable by applying a magnetic field along the anosotropy axis. 
Burning, of course, should lead to a much faster relaxation than a regular 
relaxation at fixed low temperatures. Indeed, in early experiments on relax- 
ation of large specimens of MM [5jj7] an abrupt and nearly total relaxation 
of the metastable magnetization has been detected but not explained. The 
2005 space-resolved experiments of the Sarachik group [8] on long crystals 
of Mni2Ac have shown propagating fronts of relaxation. In this experiments, 
regularly-spaced Hall probes at the sides of the crystal detected the trans- 
verse magnetic field created by the non-uniformity of the magnetization [9]. 
Chudnovsky interpreted these propagating fronts of relaxation as fronts of 
deflagration p\. Measurements of the time dependence of the total magneti- 
zation by the Tejada group, inspired by the above experiment, have shown 
a linear time dependence that was attributed to a deflagration front travel- 
ling through a Mn 12 crystal |10|. Here quantum maxima of the front speed 



vs the bias field have been detected, Fig. 4 of Ref. 10 . Discovery of mag- 
netic deflagration opened an active field of experimental research, mainly 
on Mni2Ac |Tl T4 . Experiments at high sweep rates (15) [16] have shown 
spin avalanches propagating at a fast speed. In this region, deflagration can 
go over into detonation [17]. Magnetic deflagration (coupled to a structural 
phase transition) has also been observed on manganites [18] and intermetallic 
compounds jl9][20]. To the contrast, it is problematic to observe deflagration 
fronts on another popular MM Fes because of the pyramidal shape of its 
crystals. 

One can ask if deflagration can exist in traditional magnetic systems, many 
having a strong uniaxial anisotropy. Unfortunately, the energy release in mag- 
netic systems is much weaker than in the case of a regular (chemical) defla- 
gration. Thus, at room temperatures, the ensuing temperature increase is 
too small to change the relaxation rate and support burning. Only at low 
temperatures the increase of the relaxation rate becomes large. A hallmark 
of magnetic deflagration is its non-destructive character. "Burned" MM can 
be recycled (put again into the metastable state) by simply reversing the 
longitudinal magnetic field. 



A comprehensive theory of magnetic deflagration given in Ref. 21 includes 
calculations of the stationary speed of the burning front, ignition time due to 
local increase of temperature or change of the magnetic field, as well as the 
analysis of stability of the low-temperature state with respect to deflagration 
that depends on the heat contact of the MM crystal with the environment. 
However, up to now there is no complete accordance between the theory 
and experiment for several reasons. First, thermal diffusivity k of Mni2 that 
plays a crucial role in deflagration has not been measured up to now. Second, 
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there is no completely satisfactory theory of relaxation in molecular magnets 
that takes into account important collective effects such as superradiance and 
phonon bottleneck. 

Because of their not too large spin (S — 10 for Mni2 and Fes), molecular 
magnets arc famous exponents of spin tunneling 22 - 25 that has a resonance 
character and leads to the steps in dynamic hysteresis curves at the values 
of the longitudinal magnetic field where quantum levels of the spin at the 
two sides of the potential barrier match [26 ■ 28 . Since the discovery of mag- 
netic deflagration there was a quest for quantum effects in it. The simplest 
approach (lO|[2T] uses the fact that usually spin tunneling occurs via pairs of 
quantum levels just below the classical barrier. This tunneling is thermally 
assisted and can be described by an effective lowering of the energy barrier 
at resonance values of the bias field (Fig. 2 of Ref. [7]). Thus using the Ar- 
rhcnius relaxation rate with such an effective barrier does incorporate spin 
tunneling. Experimentally it was found that spin tunneling strongly affects 
ignition of deflagration (Fig. 5 of Ref. (TTj) and to a smaller extent the front 
speed (Fig. 5 of Ref. [IT] and Fig. 4 of Ref. (To]). 

Quantum effects in deflagration shoud be sensitive to the dipolar field 
created by the sample. In a long uniformly magnetized crystal of Mn 12 Ac 
the dipolar field is B^ = 52.6 mT, as calculated microscopically in Ref. 
[29] , while the measured value [14] is very close to it. This creates a dipolar 
energy bias = g/isB^ (m' — m) between the pair of resonant quantum 

levels m and m (quantum numbers for S z in the two energy wells). This 
energy bias typically largely exceeds the tunnel splitting A that contributes 
to the resonance width. In the deflagration front the dipolar field typically 
changes between +B( D > and —B^ and so does the energy bias. As the result, 
spin tunneling in the deflagration front does not occur at a fixed resonance 
condition. This can explain why the observed quantum maxima in the front 
speed can be not as strong as expected, compared to the effect of tunneling 
on the ignition of deflagration. 

Further theoretical research led to the idea of the dipole-dipole interaction 
(DDI) playing an active role in deflagration by controlling the relaxation 
rate, as temperature does in the regular deflagration. Adding to the external 
bias field, the dipolar field can set particular magnetic molecules on or off 
resonance, facilitating or blocking their tunneling relaxation. The problem is 
self-consistent since tunneling of one magnetic molecule changes dipolar fields 
on the other ones. A numerical solution of this problem in a form of a moving 
front of tunneling at zero temperature (sometimes called "cold deflagration") 



has been found in Ref. 30 . An analytical solution for the front of tunneling 
in the realistic strong-DDI case has been obtained in Ref. [31] . 

Pure non-thermal fronts of tunneling can occur in the case of a very good 
thermal contact of the MM crystal with the environment, so that its temper- 
ature does not increase and remains so low that tunneling takes place directly 
from the metastable ground state into a matching excited state on the other 
side of the barrier. This process can be efficient only if a strong transverse 
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field is applied and the corresponding tunnel splitting A is large enough. In 
this case the speed of fronts of tunneling can theoretically exceed the speed 
of a regular deflagration by a large margin. Indeed, the dipolar field in the 
crystal changes instantaneously, in contrast to the temperature changing via 
heat conduction. Second, the relaxation rate due to tunneling directly from 
the ground state can be much higher than the relaxation rate due to the 
barrier-climbing processes in the regular deflagration. 

If the MM crystal is thermally insulated, its temperature is increasing as 
a result of a decay of the metastable state, so that there can be a mixture 
of both mechanisms of deflagration considered above [32] . Whereas far from 
resonances a regular deflagration takes place, near resonances tunneling leads 
to a great increase of the front speed. A more detailed treatment of the 
quantum-thermal deflagration for a realistic model of Mn 12 Ac with S% terms 



in the effective Hamiltonian is given in recent Ref. 33 



Theories of fronts of tunneling mentioned above are based on the model 
simplification considering it as one dimensional. In the regular deflagration, 
there is a mechanism that makes fronts flat and smooth (laminar), so that 
the deflagration problem in long crystals indeed becomes Id. In the case of 
dipolar-driven fronts of tunneling, it is not immediately clear whether fronts 
are flat or not, and, moreover, there is a mechanism that favors non-laminar 
fronts. The full 3d theory of fronts of tunneling that will be presented below, 
numerically yields non-flat and non-laminar fronts. The latter slows down 
the front speed in comparison to the simplified Id theory but, nevertheless, 
the speed can reach values comparable with the speed of sound in MM near 
tunneling resonances in strong transverse field. 

In the main part of this contribution, first the regular (thermal) mag- 
netic deflagration will be considered. Then calculation of the dipolar field in 
molecular magnets will be explained. The final part is devoted to the theory 
of fronts of tunneling. 



1.2 Magnetic deflagration 

For the generic model of a molecular magnet the energy has the form 

H = -DS 2 Z - gfi B B z S z + W, (1.1) 

where D > is the uniaxial anisotropy constant and %' stands for all terms 
that do not commute with S z and thus cause spin tunneling. In Mn^Ac there 
is an additional smaller longitudinal term ~AS Z , the implications of which 
will be discussed later. In the biased case B z > 0, the dependence of the 
energy on <j z = S z /S is sketched in Fig. [TTj The energy barrier U shown in 
Fig. 1 1 - 1 1 has the form 
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Fig. 1.1 Energy barrier of a biased molecular magnet, <r z = S z /S. 



U = (l 



h) 2 U , 



U a = DS 2 



h = gn B B z /(2DS). (1.2) 



With S = 10 the zero-field energy barrier U$ has a large value of 67 K in 
Mni 2 Ac. The energy of the metastable state is given by AE = 2Sg[iBB z . 

In the absence of spin tunneling at low temperatures, U / (ksT) 3> 1, 
the rate equation describing relaxation of the metastable population n (the 
fraction of magnetic molecules in the left well) has the form 



where the relaxation rate is given by 



r = r exp 



U 



1 + exp 



AE 

'knf 



(1.3) 



(1.4) 



Here the second term in the square brackets describes back transitions from 
the stable well to the metastable well. In the strong-bias case, AE ^> ksT, 
this term can be omitted. The equilibrium metastable population 

n (eq) is 

given by 



,(eq) _ 



= V 



exp 



AE 



1 



(1.5) 



In the strong-bias case it can be neglected. 

The second equation describing deflagration is the heat conduction equa- 
tion 

CT = V • fcVT — nAE, (1.6) 
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where k is thermal conductivity and C is heat capacity. The second term 
on the right is the energy release due to decay of the metastable state. The 
heat capacity is mainly due to phonons, whereas the magnetic contribution 
is relatively small. At low temperatures only acoustic phonons are excited, 
whereas high-energy optical phonons are frozen out, thus C has the form |34| 

C = Ak B (T/0 D ) a , (1.7) 

where a — 3 in three dimensions, A = 127r 4 /5 ~ 234 is a numerical factor 
and &£> is the Debye temperature, Op ~ 40 K for Mn 12 Ac. Although at low 
temperatures this expression is in a reasonable accordance with measure- 
ments on Mni2Ac [35 , its applicability range is very narrow, T < 5 K. On 
the other hand, the temperature generated in the deflagration (the so-called 
flame temperature) is typically above 10 K. The heat capacity of Mn^Ac 
can be well described within a broad temperature range with the help of the 
extended Debye model (EDM) |36| that comprises three different acoustic 
phonon modes as well as optical modes. Practically, one can use measured 
values of C 135 . 



It is convenient to use the relation C — d£/dT to rewrite Eq. (1.6) in 
terms of the energy £ as 

£ = V • nV£ - nAE, (1.8) 

where k — k/C is thermal diffusivity. The latter has not yet been measured, 
although a crude estimate k ~ 10~ 5 m 2 /s was deduced from experiments 



8|13 . This value is comparable with that of metals. Temperature dependence 



of k that could be substantial at low temperatures remains unknown. 



Equations (1.3) and (1.8 1, together with Eq. (1.4 1 and the relation 



£{T) = / C(T')dT', (1.9) 
Jo 

is a strongly-nonlinear system of equations. It is easy to solve these equations 
numerically but it costs efforts to do it analytically. The two main problems 
to solve are (i) stability of the low-temperature state with respect to thermal 
runaway or ignition of a deflagration front and (ii) the shape and speed of 
the stationary deflagration front in long crystal. 



1.2.1 Ignition of deflagration 

If the sample is perfectly thermally insulated, the whole released energy re- 
mains inside and the temperature monotonically increases. This leads to a 
thermal instability that can take a considerable time to develop, the ignition 
time. If there is a thermal contact with the environment, maintained at a 
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Fig. 1.2 Semenov's mechanism of a thermal runaway, Eq. ( 1.10 i 



constant low temperature To, there are two possible cases. In the subcritical 
case, the temperature rise in the sample due to slow decay leads to temper- 
ature gradients and heat flow toward the sample boundaries that ensures a 
stationary low-temperature state (proper conditions of explosives' storage). 
In the supercritical case, heat loss through the boundaries is insufficient to 
balance the increase of the heat release due to rise of temperature. This 
leads to ignition of a self-supporting burning process. In small crystals of 
MM, temperature gradients are higher and heat loss to the environment is 
more efficient. In larger crystals, temperature gradients are lower and ther- 
mal instability is more likely. This is why deflagration was observed in larger 
crystals. 

Thermal instability occurs because of a stronger temperature dependence 
of the relaxation rate, Eq. (1.4 1, than that of the heat exchange with the 



environment. The essence of the problem is contained in the old model of 
explosive instability by Semenov described by a single equation 



T = Q r 



'cooling i 



(1.10) 



where Q rca ction ~ r(T) and Qcooiing = a(T - T ). In the case B in Fig. \l.2 
the thermal contact to the bath is too weak, Qcooiing < Qreaction at all T, so 
that the system is absolutely unstable. In the case A, the thermal contact 
is stronger and there is a stability range T < T 2 , where the stationary state 
T = T\ is an attractor. However, heating the system above T 2 leads to thermal 
explosion. 

Semenov's model is zero-dimensional, whereas in MM crystals the problem 
is at least one-dimensional and more complicated. There are different cases 
of thermal instability, mainly instability of a large crystal initially at uniform 
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temperature (that begins at the center) , instability due to heating one end of 
a long crystal, and the instability due a magnetic field gradient that makes 
the barrier lower at one side of the crystal. Analysis of all these cases has 
been done in Ref. [21] . In particular, when the magnetic field and temperature 
of the sample boundary To are independent of coordinates, the crystal loses 
stability against formation and propagation of the flame (magnetic avalanche) 
when the rate of the spin flip for an individual molecule, T(7J, To) , exceeds 

_ k B T 8kT a 
c ~ U PriiAE ' [ ' 

Here k is thermal conductivity at To and the length parameter I is uniquely 
determined by geometry, being of the order of the smallest dimension of the 
crystal, whereas m is the metastable population in the initial state. 

Experimentally magnetic deflagration can be initiated either by heating 
one end of the crystal [IX- 13 or by sweeping the magnetic field in the positive 



direction, that reduces the energy barrier and makes the condition in Eq. 
( |1.11[ ) satisfied |8|. In Ref. [lO] deflagration was ignited by surface acoustic 
waves (SAW), instead of heating. 



1.2.2 Deflagration fronts 

Fronts of magnetic burning propagating in long crystals of molecular magnets 
are flat and smooth, i.e., the problem of deflagration is one-dimensional. The 
stability of flat fronts can be immediately seen. Indeed, if a fraction of a 
front gets ahead of neighboring fractions, the heat released at this place 
will be propagating not exactly straight ahead (as in a flat front) but also 
sideways. This will slow down this leading fraction of the front and speed 
up the lagging fractions surrounding it. Thus any local deviation from a flat 
front will disappear with time. 

In a stationary-moving front, all physical quantities depend only on the 
combined variable that can be chosen, e.g., in the time-like form u = t — z/v, 
where v is the front speed. In terms of u the deflagration equations have the 
form 



-r(T) (n-n (eCL >(T)) 



dn 
du 

d£ 1 d d£ dn 
du v 2 du du du 



2 ^-.--AE (1.12) 



plus Eq. (1.9). Integrating the energy equation one obtains 



k d£ 

E + nAE ~-r= const. (1.13) 

v z du 
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Far before and far behind the front, the term with the derivative vanishes. 
Thus one obtains the energy conservation law in the form 



Si + n,AE = S f + n icq) (T f )AE 



(1.14) 



where i stands for "initial" (before the front) and / stands for "final" or 
"flame". This is a transcedental equation for the flame temperature Tf that 
has to be solved together with Eq. (1.9). If Si ~ (low initial temperature) 



and rS cq ^(Tf) is negligible (full-burning case realised at a strong bias, see Eq 
(74) of Ref. pll) one immediately finds the flame energy from riiAE = S 



and then Tf follows by inverting Eq. ( 1.9 I. In the incomplete-burning regime 



at small bias, a pulsating instability of stationary deflagration fronts |37| was 
found. The operations above assume that the heat is not exchanged via the 
sides of the crystal. In the opposite case, the energy conservation becomes 
invalid and the theory has to be extended. 

One can immediately get an idea of the front speed by rewriting the defla- 



gration equations (1.121 in the dimensionless form. In terms of the reduced 
variables 

h = n/rii, S = S/(riiAE), u = uJy (1.15) 

and parameters 

r = r/r f , k = k/kj (i.ig) 



equations (1.121 become 



du \ J 

dS 1 d _ dS dn 
du v 2 du du du 

where the reduced front speed v is related to the actual front speed v by 

v = iy/n f r } . (1.18) 

Refs. 1 2, 8 give the expression above without v for the front speed. 

It turns out that v in Eq. (1.181 is not merely a number but rather it is a 



function of dimensionless parameters such as 

W f ee U/{k B T f ). 



(1.19) 



Because of the non-linearity of Eq. (1.17), their general analytical solution 



that defines v does not exist. There are two parameter ranges in the problem: 
Slow-burning high-barrier range Wf S> 1 and fast-burning low-barrier range 
W f <l. 

In the former, burning occurs in the front region where the temperature 
is already close to Ty. Assuming that k is temperature independent, R = 1, 
and linearizing r(T) near Ty, one can solve the problem analytically. Within 
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the full-burning approximation (n( cq ) => 0) the reduced front speed is given 



by 21 



CfT f k B T f 
jiiAE U 



(1.20) 



With the help of Eq. ( 1 . 7 1 (that is not accurate, however ! ) this result simplifies 
to 

v = yJ(a+l)/W f . (1.21) 

The applicability range of these expressions is v <C 1. 

The corresponding profile of the metastable population n in the front has 
the form 

n= = 1 |l-tanh~] (1.22) 

1 + e" 2 V 2 J y ' 

that corresponds to the symmetric tanh magnetization profile a z = 1 — In = 
tanh(u/2). In real units the result reads 



1 + tanh 



2vL 



(1.23) 



where Id — y/ Kf / Tf is the a-priori with of the deflagration front. Magneti- 
zation profile of this kind can be seen in Fig. 11 of Ref. (2ll and in the upper 
panel of Fig. 10 of Ref. (33] . The reduced energy in the front is given by 



e = (i 



(i-hr 



(1.24) 



Since in the high-barrier approximation v -C 1, the formula above yields 
S 1 in the active burning region and actually everywhere except for the 
region far ahead of the front where h is very close to 1. This justifies the 
approximation made. 

It should be noted that the full-burning approximation used above requires 
a bias high enough thus the barrier low enough, Wf < 6, according to Eq. 
(79) of Ref. 21 . Thus the applicability range of the slow-burning high-barrier 



approximation is rather limited. The theory can be improved by taking into 
account incomplete burning. However, this makes analytics cumbersome be- 
cause of the transcedental equation ( 1.13 1 defining TV . Numerical solution for 



the deflagration front poses no problems, nevertheless. Because of incomplete 
burning, Tf and thus the front speed decrease below the values given above. 

In the low-barrier fast-burning regime Wf < 1 there is no rigorous analyti- 
cal solution to the problem. Additionally, the Arrhenius form of the relaxation 



rate, Eq. (1.4 1, becomes invalid. In this regime the magnetization profile is 



asymmetric, as can be seen in the upper panel of Fig. 12 of Ref. 33 
Making the symplifying approximation for the relaxation rate 



f(£) 



0, £ <£ Q 

1, £ > So, 



(1.25) 
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Reduced speed of the deflagration front 




Analytical 



10 



15 



W 20 



Fig. 1.3 Reduced speed of a deflagration front defined by Eq. 1 1.18 1. The numerical 
result has been obtained in Ref. [21j within the full-burning approximation using the low- 
temperature form of the heat capacity, Eq. ( |1.7| . Analytical result is Eq. jl.34[ | . The dotted 
line is the high-W^ asymptote, Eq. jl.21| . 



where £q will be defined below, one can solve the problem of a stationary 
deflagration front in the whole parameter range. Let us search for th e fron t 
in which £ = £q at u = 0. In the reduced form of the energy equation (1.13 1, 

^ =v 2 {£ +n-l), (1.26) 
du 

one has n — 1 before the front, u < 0. Thus here the energy equation solves 
to 

f = £ e" 2 ". (1.27) 

On the other hand, for u > the solution of the po pulati on equation dh/du = 
—Tn = —n reads n = e~ u . Inserting this into Eq. (1.261, and integrating the 
differential equation, one obtains the solution 

£=(£ --^)e^ + l--^e-\ (1.28) 
\ 1 + try 1 + v 2 - 

The first term of this expression must vanish because of the boundary con- 
dition £ (oo) = 1. This defines the reduced front speed, 



6-^-1. (1.29) 
To define £q, consider the reduced Arrhenius relaxation rate 
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r = cxp 



sr ; 1 1 - i 

T 



and require 



n , (i-il = -i 



as the switching point between t = and T = 1. This yields 



T n = 



1 + W) 



Using Eq. (1.7 1, one obtains 



Wt 



1 + W f 



a+1 



Substituting this into Eq. (1.291, one finally obtains 



, w f 



Limiting cases of this formula are 



' f 



i)/w f , w f > i 

W f < 1. 



(1.30) 



(1.31) 



(1.32) 



(1.33) 



(1.34) 



(1.35) 



It is remarkable that the rigorously obtained high-barrier slow-burning result 



of Eq. (1.21) is recovered exactly. In the low-barrier fast-burning case the 



reduced front speed becomes large, as well as the actual front speed of Eq. 



(1.18 1. One can see that Eq. (1.34) is in a good accordance with the numerical 



solution shown in Fig. |1.3| 

The high-speed regime of the deflagration should be superceded by deto- 
nation when the front speed approaches the speed of sound. In detonation, 
thermal expansion resulting from burning sends a shock wave into the cold 
region before the front. As a consequence, the temperature before the front 
rises as a result of compression, initiating burning. Such a mechanism was 
recently considered for Mn^Ac in Ref. [17) . 



1 Theory of deflagration and fronts of tunneling in molecular magnets 



13 



1.3 Fronts of tunneling 

1.3.1 Tunneling effects in the relaxation rate 

The relaxation rate r including spin tunneling is at the foundation of the 
quantum theory of deflagration in molecular magnets. In the generic model 



of MM, Eq. ( 1.1 ), tunneling resonances occur at the values of the total bias 



field B tot z (including the self-produced dipolar field) equal to 

B k = kD/(gvB), fc = 0,±l,±2,... (1.36) 

for all the resonances. Spin tunneling leads to the famous steps in the dynamic 
hysteresis curves [ 26-28 . In the real Mn 12 Ac there is an additional term 



-ASl; that makes higher-energy resonances be achieved at smaller B z than 
low-energy resonances. The resulting tunneling multiplets 

9H B B km = k[D+ (to 2 + (to + kf) A] (1.37) 



were used to experimentally monitor 38|39 the transition between thermally 



assisted and ground-state tunneling 40 in Mn 12 Ac. Below B k will stand for 



the resonance field B km , for simplicity of notations. 

In the case of an isolated magnetic molecule, the probability of a spin to 
be in one of the resonant quantum states is oscillating with time with the 
frequency A/h, where A is the tunnel splitting. However, coupling to the 
environment, e.g., to phonons, introduces damping to these oscillatins. If the 
decay rate of at least one of the resonance states, F m or r m >, exceeds A/h, 
tunneling oscillations of the spin are overdamped. This can be illustrated in 
the case of a resonance between the metastable ground state \—S) and the 
matching excited state at the other side of the barrier \m!) of a biased MM 
at zero temperature. Ignoring all other levels, that is justified at T — 0, one 
can write down the Schrodinger equation in the form |31| 

* A 

W l r \ i A 

-Y~2 rm ') Cm '~2h C ~ S ' (08) 
where 

W = e- s - em' =(S+ m')gnB(B t ot,z - B k ) (1.39) 

is the energy bias between the two levels. Whereas the level |— S) is un- 
damped, the level \m') can decay into lower-lying levels in the same well 
via phonon-emission processes. At T — there are no incoming relaxation 
processes for \m'). In this case the damped Schrodinger equation above is 
accurate, as it can be shown to follow from the density matrix equation. In 
the underdamped case r m i < A/h the solution of these equations is oscillat- 
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ing. The first choice for studying tunneling dynamics in molecular magnets 
is the overdamped case r m > ^> A/h, since for not too strong transverse fields 
B± the tunnel splitting/! is a high power of B± (Ref. 141] ) and typically it 
is much smaller than r m >. In the overdamped case the variable c m ' in Eq. 



(1.381 adiabatically adjusts to the instantaneous value of c_s and the solu- 
tion greatly simplifies. Setting c m < = in the second of these equations, one 
obtains 

Cm ' = 2hw/h + ir m ,/2- (L40) 



Inserting this into the first of equations (jl .38|) yields a closed differential 
equation for c_g. Using n = | cl- 
one arrives at the rate equation 



equation for c_g. Using n — |c_g| 2 for the metastable occupation number 



-rn, (1.41) 



where the dissipative resonance-tunneling rate r is given by 42 



a 2 r ill 

r =^7 2- (!- 42 ) 

2^ 2 {w/hf + (r m/ /2) 2 

This is a Lorentzian function with the maximum at the resonance, W = 0. 



Eqs. (1.411 and (1.42) were used in Refs. 30 31 to study dipolar-controlled 



fronts of tunneling at T = 0, or "cold deflagration". The full system of equa- 



tions ( 1.38 ) could also be used to this purpose but nothing had been published 
up to date. 

At nonzero temperatures, tunneling transitions via higher energy level 
pairs become possible (thermally-assisted tunneling) and one has to take 
into account non-resonant thermal transitions over the top of the barrier. 
This makes the problem more complicated, and one needs to use the density 
matrix equation (DME) taking into account spin-phonon interactions explic- 
itly. One of the first works using DME for Mn 12 Ac was Ref. |42 in which 
spin tunneling was considered with the help of the high-order perturbation 
theory [IT] for a small transverse field B± . The spin-phonon processes taken 
into account were due to dynamic tilting of the anisotropy axis by transverse 
phonons. Ref. |42| could qualitatively explain thermally-assisted tunneling 
via the level pairs just below the classical barrier. However, tunneling via 
low-lying resonant level pairs or tunneling directly out of the metastable 
ground state are inaccessible by this method because large enough splitting 
requires non-perturbatively large transverse field that can only be dealt with 
numerically. 

Further work on spin-phonon relaxation in MM lead to elucidation of the 
universal relaxation mechanism |43[|44| . This mechanism consists in distor- 
tionless rotation of the crystal field acting on a magnetic molecule, actually 
the same mechanism as used in Ref. (421 . It was, however, understood that this 
mechanism does not require any poorly-known spin-lattice coupling constants 
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and everything can be expressed through much easier accessible crystal-field 
parameters. This mechanism was overlooked in older theories of spin-lattice 
relaxation based on energy contributions responsible for the coupling. Ro- 
tations, to the contrary, cost no energy and the effect has a purely inertial 
origin. 




3 B z (T) 



Fig. 1.4 Relaxation rate of Mn^Ac vs temperature and longitudinal magnetic field in a 
small transverse field. Resonance multiplets with k = 2, 3 are seen. 




Fig. 1.5 Relaxation rate of Mni2Ac vs temperature and longitudinal magnetic field in 
the transverse field B± = 3.5 T. One can see the ground-state resonance at B z = 0.522T 
and the first-excited-state resonance at B z = 0.490T for k = 1 multiplet. 
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The universal relaxation mechanism allows a general numerical implemen- 
tation of the DME fully based on the crystal field parameters, recently ac- 
complished in Ref. 45 that summarizes the current state of the problem. 



Another important feature of Ref. 45 is using the so-called semi-secular ap- 



proach capable of dealing with resonant pairs of levels and thus describe spin 
tunneling. Conventional implementations of the DME (see, e.g., Ref. [46]) use 
the secular approach that crashes on tunneling resonances. In Ref. |45| the 
relaxation rate r is extracted from the time-dependent numerical solution of 
the DME (expressed in terms of eigenvalues and eigenfunctions of the density 
matrix) as the inverse of the integral relaxation time (47p8 . Unlike using the 
lowest eigenvalue of the density matrix, this method also works at elevated 
temperatures. 

The temperature and field dependence of r in Mn 12 Ac at a small trans- 
verse field (B± = 0.04T that typically arises due to a 1 Q misalignment of 
the easy axis and the applied longitudinal field) is shown in Fig. 1.4 One 
can see very narrow and high maxima of _T (note that log r is plotted!) due 
to spin tunneling. Maxima corresponding to the ground-state tunneling, for 
which the maximum in r does not disappear at T = 0, correspond to the 
highest value of B z in the multiplet. There are k = 2 and k = 3 tunneling 
multiplets seen in Fig. |1.4| Note that tunneling via low-lying resonances is 
relatively weak and it is eclipsed by the thermal activation contribution at 
higher temperatures. 

At stronger transverse field such as B± = 3.5T in Fig. |1.5| the barrier is 
strongly reduced and high-lying tunneling resonances are broadened away. 
Here, one can see the ground-state resonance (B z = 0.522T) and the first- 
excited-state resonance (B z = 0.490T) for k — 1 multiplet. The ground- 
state resonance does not disappear at the highest temperature that has an 
important implication in the dynamics of fronts of tunneling. Note the much 
higher tunneling rate at T — 0, in comparison to the previous figure. 

A puzzle in the theory of relaxation of molecular magnets is the prefac- 
tor Jq in the Arrhenius relaxation rate, Eq. (1.4 1, being by two orders of 



magnitude too small. This was already recognized in the early work 42 



Using the standard spin-lattice relaxation model considering one spin in an 



infinite elastic matrix, it is impossible to arrive at r 
in experiments 



10 s observed 



interactions 



50 



35 49 



without introducing artificially strong spin-phonon 
For a strongly diluted molecular magnet, considering a sin- 
gle spin could be justified, but in the regular case it can not. High density 
of magnetic molecules should lead to such collective effects as superradi- 

Possibility of superradiance in 



51-53 and phonon bottleneck 54 -56 



ance 

fast avalanches triggered by a fast field sweep has been discussed in Ref. [15] 
References 57, 58 report microwave emission from MM that can be inter- 
preted as superradiance. However, it would be difficult to address these com- 
plicated issues while dealing with the quantum deflagration problem, so that 
the calculated relaxation rate will be simply multiplied by 100 to approxi- 
mately match the experiment, as was done in Ref. 1 33 . 
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I. 3.2 Dipolar field in molecular magnets 

Very sharp resonance peaks in the relaxation rate r seen in Figs. |1.4| and 

I I. 5 1 require an accurate calculation of the dipolar field in the crystal that can 
self-consistently control tunneling by setting individual molecules on or off 
resonance. The equations describing this are the same relaxational equation 



(1.3 1 and thermal equation (1.8), as before, only with r depending on the 
total magnetic field 

B, ot|J (r)=B ? + Bf(r), (1.43) 

where B z is the external bias field and B Z D ^ is the self-consistently calculated 
dipolar field. In the case of cold deflagration, the thermal equation can be 



discarded and one has to solve only the relaxational equation (1.41 1. Since 
the dipolar field depends on the magnetization everywhere in the crystal, the 
equations of quantum deflagration are integro-differential equations. Note 
that the transverse component of the dipolar field can be discarded because 
its effect is small. 

For the purpose of calculating the dipolar field, conventional magneto- 
statics (see, e.g., Ref. [59]) is unsuitable because it provides an irrelevant 
magnetic field formally averaged over the microscopic scale that ignores the 
lattice structure. The physically relevant dipolar field is the field created at 
positions of magnetic molecules by all other molecules. It is a microscopic 
quantity that depends on the lattice structure. To illustrate this point, mag- 
netostatic field in a uniformly magnetized long sample is — AttM. where 
M is the magnetization. However, microscopically calculated dipolar field in a 
long uniformly magnetized crystal of Mn^Ac is much smaller, B Z D ^ = 5.26M. 

It is convenient to express the z component of dipolar field at site i (i.e., 
at a particular magnetic molecule) in the form 

BW = (Sgii B /vo)D„, (1.44) 

where vq is the unit-cell volume. For Mn 12 Ac one has SgfiB/ v o — 5.0 mT. 
The reduced dipolar field D zz , created by all other molecular spins polarized 
along the z axis is given by 

r> V^j. j 3(e z • n 4J ) 2 - 1 r i:j 

Di,zz = Z^VijVjz, Pij = vo 3^ , n tJ = (1.45) 

where a z = S z /S. To calculate the sum over the lattice for the site i, one 
can introduce a small sphere of radius vq around i satisfying v^ 3 < ro < I, 
where L is the (macrocopic) size of the sample. The field from the spins at 
sites j inside this sphere can be calculated by direct summation over the 
lattice, whereas the field from the spins outside the sphere can be obtained 
by integration. The sum of the two contributions does not depend of tq. If 
the magnetization in the crystal depends only on the coordinate z along the 
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symmetry axis of the crystal that coincides with the magnetic easy axis z 
(that is the case for a fiat deflagration front), the integral over the volume 
can be expressed via the integral over the crystal surfaces. The corresponding 
contribution can be interpreted as that of molecular currents flowing on the 
surface. The details are given in the Appendix to Ref. 



29 



In particular, for a uniformly magnetized ellipsoid the total result has the 
form 



3 



kj = D zz a z , 



independently of i, where 



D zz = £>^ ph) + Airv fl/3 - t»W) 



(1.46) 



(1.47) 



and v is the number of molecules per unit cell (v = 2 for Mn^Ac having 
a body-centered tetragonal lattice). Here comes from the summation 

over a small sphere and the remaining terms come from the integration. For 
the demagnetizing factor one has rS z ^ — 0, 1/3, and 1 for a cylinder, sphere, 



and disc, respectively. One obtains D 



(sph) 



D 



(sph) 



< for a tetragonal lattice with a 



= for a simple cubic lattice, 
b < c, and L»iT h) > for that 



with a = b > c. The latter is the case for Mn^Ac having D 
a long cylinder this results in £j|z = 10.53 or, in real units 

B{ D} = 52.6 mT. 



(sph) 



2.155. For 



14 



2!) 



(1.48) 



The dipolar energy per magnetic molecule can be written in the form 
E = -(1/2)D ZZ E D , where 



E 



D 



(Sgfj, B ) A'o 



(1.49) 



is the characteristic dipolar energy, Ejy/ks = 0.0671 K for Mn^Ac. The role 
of the DDI in spin tunneling is defined by the ratio of the typical dipolar bias 



W<- D) = 2Sg/i B Bi D) = 2E D D y zz yl > to the width of the overdamped tunneling 
resonance r m i in Eq. ( 1.42 I. It is thus convenient to introduce the parameter 



,(cyl) 



E D = 2E D /(hr m ,) 



(1.50) 



that is always large. For instance, using the experimental Arrhenius prefactor 
r ~ 10 7 s _1 for r m t. one obtains E D ~ 10 3 . 

For a cylinder of length L and radius R with the symmetry axis z along 
the easy axis, magnetized with a z = & z (z), the reduced dipolar field along 



the symmetry axis has the form 29 



D zz {z) 



L/2 ^ 2^R 2 a z (z') 
L ' 2 \{z'-z) 2 + R 2 



3/2 



k D o z {z), 



(1.51) 
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Fig. 1.6 Reduced dipolar field in a deflagration front in the slow-burning limit, created 
by the magnetization profile <r z (z) = — tanh [(z — zq)/Ij\. Analytical result: Eq. | |1,51| ; 
Points: Direct summation of dipolar fields over Mni2Ac lattice. 



where a z — 1 — 2n is polarization of pseudospins representing spins of mag- 
netic molecules (a z = ±1 in the ground and metastable states, respectively) 
and 

(1.52) 



k D = 87TZ//3 - L>i s z ph ) = Attv - D<g® > 0, 



fcrj = 14.6 for Mn 12 Ac. In Eq. (1.511, the integral term is the contribution of 
the crystal surfaces, while the lattice-dependent local term is the contribu- 
tion obtained by direct summation over lattice site within the small sphere 
ro minus the integral over this sphere that must be subtracted from the in- 
tegral over the whole crystal's volume. For other shapes such as elongated 



rectangular, one obtains qualitatively similar expressions 31 



A striking feature of Eq. (1.511 is that the integral and local terms have 
different signs. The integral term changes at the scale of R while the local term 
can change faster, that creates a non-monotonic dependence of D zz (z). In the 
case of a regular magnetic deflagration, the spatial magnetization profile in 
the slow-burning limit is of the type a z {z) = — tanh [{z — zq)/^], where Id 
is the width of the deflagration front that satisfies <C R, c.f. Eq. (1.231. 
The resulting dipolar field is shown in Fig. |1.6[ where the line is the result of 



Eq. ( 1.51 1 and points represent the dipolar field along the symmetry axis of a 



long cylindrical crystal calculated by direct summation of microscopic dipolar 



fields over the Mni^Ac lattice. One can see that Eq. (1.51 1 is pretty accurate, 
small discrepancies resulting from being not large enough in comparizon 
to the lattice spacing a. The central region with the large positive slope 



is dominated by the local term of Eq. (1.51) that changes in the direction 
opposite to that of the magnetization. For R ^> Id, D zz reaches the values 
±14.6 due to the local term before it begins to slowly change in the opposite 
direction. In real units the dipolar field at the local maximum and minimum 
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is ±B z k °\ where 

B (k D ) _ 72.9 m T, (1.53) 
exceeding the dipolar field of the uniformly magnetized long cylinder Eq. 



( 1.48 1. Also one can see from Fig. 1.51 that the dipolar field becomes opposite 



to the magnetization at the ends of the cylinder, that should lead to an 
instability of the uniformly-magnetized state in zero external field. 

The Id theory of fronts of tunneling 30 33 is based on the simplifying 



assumption that the deflagration front is flat, a z = a z [z), and the dipolar 



field is given by Eq. (1.51) everywhere. Since, in fact, the dipolar field also 
depends on the distance from the crystal's symmetry axis, it is likely that 
such a more complicated structure of B z will self-consistently affect the front 
structure, making it non-flat. 

There is also a question of stability of a smooth front at a small scale. 
Whereas flat and smooth fronts of regular burning are stable, there is an 
instability mechanism for a flat front in the presence of tunneling controlled 
by dipolar fields that will be explained below. This is why it is important to 
develop a full 3d theory of fronts of tunneling. 

If the magnetization a z of a MM crystal depends on all the coordinates 
x, y, z but this dependence still has a macroscopic scale, one can again use the 
method of calculating the dipolar field that combines summation over a small 
sphere (where a z does not change) and integration over the remaining volume 
of the crystal. In this case the integral over the volume does not reduce to an 
integral over the surface and it has to be done numerically. In the solution 
of the deflagration problem, it is convenient to discretize the volume of the 
crystal and use the same grid to sample the magnetization variables and to 
calculate the dipolar field. A problem with this integral is that a small sphere 
of radius tq (around earch point r = r,; where the dipolar field is calculated) 
has to be excluded from integration and the contribution of this excluded 
region is comparable with the total result because of the singularity of the 
DDI. 

The solution to this problem is, for any point r, to add and subtract the 
dipolar field in a uniformly magnetized crystal with <j z = er z (r). The total 
reduced dipolar field can be thus represented as 

D zz {r)= j dv'4>{v' -v){a z {v r )-a z {v))+a z {r)(V zz {v)-k D ), (1.54) 



where <f> is defined in Eq. ( 1.45 1. Because of the terms subtraction at r' — > r, 
the contribution of the excluded small sphere in the intcrgal is negligible and 
the integral can be extended to the whole volume of the crystal. Then the 
values of the integral for all points of a rectangular grid can be computed via 
a summation method based on the fast Fourier transform (FFT) that takes ~ 
N\og(N) operations, where N is the number of grid points. Straightforward 
calculation of the integral costs ~ N 2 operations and it has to be avoided. 
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The remainder of Eq. ( 1 .54 1 corresponds to a uniformly magnetized crystal 



and its structure is similar to Eq. (1.51 1. Again, the term with fcrj is the local 
contribution, while T> zz (r) is the contribution of surface molecular currents, 
the result of conventional magnetostatics. For a crystal of a rectangular shape 
with dimensions L x X Ly X L z the result can be obtained as a particular case 
of Eq. (88) of Ref. |60] and it has the form 



f zz( r ) = arctan 



(L x + r) x x) (L y + r] y y) (L z + r] z z) 
\J (Lx + Vxx) 2 + (L y + r] y y) 2 + (L z + r\ z z) 



-+(x = 
(1.55) 



in total 16 different arctan terms. 



1.3.3 Fronts of tunneling at T = 



The theory of dipolar-controlled fronts of tunneling at T = ("cold defla- 
gration") (30j[31 uses the relaxational equation (jl .411) with the resonance 
tunneling rate of Eq. (1.421, in which the energy bias W is given by Eq. 
(jl .39 1 with B to t,z of Eq. (1.431. Within the Id approximation 



30 



31 



the 



dipolar field is given by Eqs. (|1.44| and (|1.51|) for a cylinder. The problem is 

W/ (2E D ) that has 



thus an integro-differential equation. 

It is convenient to use the reduced energy bias W 
the form 



n- = ir.. xt +L> : ... \v^ = {S+ li )9l * B (B z -B k ), 



2E 



D 



(1.56) 



where m' = S — k is close to S for not too strong bias. Propagating dipolar- 
controlled fronts of tunneling have been found numerically |30[|31| and ana- 



lytically 31 within the dipolar window near the resonance 

0< Wext <Dt l \ 



(1.57) 

where D z c z y ^ = 10.53. In real units this yields the dipolar window 

B k < B z < B k + Bi D \ (1.58) 



where Bz is given by Eq. (|l.48|) for Mn 12 Ac 



The solution for the front of tunneling depends on several parameters such 
as the trans verse size of the crystal R and the resonant value of the relaxation 
rate of Eq. (1.421, r ics = A 2 /(h 2 r m i). Rewriting the equations in a reduced 
form j3l], one immediately finds that t he fro nt speed is of order r res R. The 
only non-trivial parameter is Ed, Eq. (1.501. An analytical solution of the 



problem is possible because of the large value of Ed. The front speed is given 
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by |3 

* ' 

v = v*T ICS R, v* ~ z —r^- ) , (1.59) 

Bk + B z — B z 

within the dipolar window, independently of Ejj. Above Bk + B Z D ^ the front 
speed is zero. The reason for this is that for the external field above Bk + 
B Z D ^ , the total field well before the front (where all spins are directed in the 
metastable negative direction and produce the dipolar field — B Z D ^) is above 
its resonance value Bk (and spin tunneling would even increase the total 
field). Thus in this case resonance tunneling transitions cannot occur. To the 
contrast, just below Bk + B Z D ^ the field well before the front is a little bit 
below the resonance and increases closer to the front where the magnetization 
is switching. In this case, there is a wide region where the system is close to 
the resonance, and the front speed becomes very high. Thus as B z crosses 
the value Bk + B Z D ^ from below, the front speed diverges and then drops 
abruptly. 

Let us compare the speed of fronts of tunneling v ~ r ics R with the speed 



of regular deflagration, Eq. (1.181. With a sufficiently strong transverse field 
applied, one can have A/h ~ r m i at the applicability limit of the overdamped 
approximation, and then r xes ~ r m i ^> J 1 / because thermal activation goes 
over high levels of the magnetic molecule where the distances between the 
levels and thus the energies of phonons involved are much smaller than for 
the low-lying levels, and also because i~y is exponentially small since Tf < U. 
Additionally, estimation of Z<j with Kf = 10 -5 m 2 /s and the experimental 
value r = 10 7 s -1 yield l d ~ 3 x 1(T 4 mm for B z near the first tunneling 
resonance and even smaller for larger bias. As in the experiment the width of 
the crystal was much larger than Id (0.3 mm in Ref. [8], 0.2 mm in Ref. 11 



and 1 mm in Ref. |10| ), one can see that r ies R 3> -T/Zd is quite possible 
in a strong transverse field, and then the front of spin tunneling is much 
faster than the front of spin burning. A very conservative estimation with 
r Tes =S> r = 10 7 s _1 and v* =*> 1 for the crystal 0.2 mm thick yields v — 1000 
m/s. As said above, in a strong transverse field one can have _T rcs 3> /bj so 
that the speed of a spin-tunneling front can easily surpass the speed of sound 
that is about 2000 m/s in molecular magnets (see analysis in Ref. |36 ). 

A hallmark of the cold deflagration is residual metastable population be- 
hind the front |31| that can be rewritten as 

n f = (B z -B k )/BW (1.60) 

(here n = nt = 1 before the front). One can see that the change of n across 
the front An = 1 — nt goes to zero at the right border of the dipolar window, 

B z = Bk + B Z D ^ . This reconciles the situation with the general requirement 
that the rate of change of the magnetization of the crystal M, limited by the 
tunneling parameter A, remains finite. Indeed, 
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Fig. 1.7 Spatial profiles of the metastable population n and the reduced bias W in the 
front for Wcxt = 2 and Ejy = 20. Everywhere in the front the system is near the resonance, 
W 0. At this value of Wext the solution begins to lose stability and periodic structures 
behind the front begin to emerge. 

M oc (1 - n f )v = r rcs R(B z - B k ) /B^ (1.61) 

reaches only a finite value M oc r ics R at the right border of the dipolar 
window before it drops to zero. 
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Fig. 1.8 Reduced front speed v* of Eq. ( 1 .59 i vs the reduced bias Wcxt of Eq. (1.56 I for 



different number of grid points. For Wcxt ^ 1 (the laminar regime) the numerical results 
are in a good accordance with Eq. l|1.59[ (straight line). 
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To obtain a numerical solution for the cold deflagration, the integro- 
differential equation was discretized to make the integral in Eq. (1.51) a 
sum and the whole problem a set of coupled non-linear first-order differential 
equations. The program was written in Wolfram Mathematica. A typical re- 
sult for spatial profiles of the metastable population n and and total energy 



bias W are shown in Fig. 1.7 In the cold deflagration front, magnetization 
and dipolar field are self-consistently adjusting in such a way that inside the 
front core of the width R the spins are on resonance and can tunnel. To the 
contrary, before and after the front magnetic molecules are off resonance and 
tunneling is blocked. One of the reasons why fronts of tunneling can be so 
fast is that their width R entering the expression for the front speed, Eq. 
(1.591, is much is much larger than the width of the deflagration front Id, 



c.f. Eq. (1.181. The solution shown in Fig. 1.8 is an example of the laminar 



solution for the cold deflagration front that is realized for a not too strong 
bias, W ext < 1-2 or B x - B k < 5-10 mT. 

For a stronger bias, the laminar solution becomes unstable. The front of 
tunneling is moving with a non-constant speed, leaving spatially-nonuniform 
distribution of the unburned metastable population behind. The spatial de- 
pendence of the dipolar field becomes discontinuous and the resonance con- 
dition in the front is not fulfilled (see Fig. 6 of Ref. (3l|). As a result, the 
front speed begins to decrease as the instability develops with the increase of 



the bias, Fig. 1.8 The instability of the solution is manifesting itself in the 
dependence on the discretization, absent in the laminar regime. 

The only experimentally feasible method to ignite cold deflagration is the 
sweep of the bias field B z . When B z is swept in the positive direction in a 
negatively magnetized MM crystal, the resonance condition is first achieved 
at the ends of the crystal where the (negative) dipolar field is weaker (see, 



e.g., the right side of Fig. 1.6). Spin tunneling at the ends of the crystal 



caused by field sweep leads to change of the dipolar field in this region that 
brings the system closer to the resonance in a region of the depth of order 
R, the transverse size of the crystal. At some moment, a spatial structure 
close to a stationary front of tunneling is formed and it begins to propagate 
into the depth of the crystal, the field sweep playing no role anymore. This 



mechanism is illustrated in Fig. 9 of Ref. 31 . Numerical calculations show 
that front of tunneling is ignited at the "magic" value of the reduced bias 
Wext — 5, weakly dependent on Ee>. For this value of the bias, the front of 
tunneling is non-laminar. 

Fronts propagating at other values of the bias, including laminar fronts, 
can be ignited by a modified procedure proposed in Ref. |31| . First, a global 
bias is being changed, as before, by a uniform field sweep until the desired 
value of W ex t is reached. After that, front of tunneling can be ignited by a 
local increase of the bias near the crystal's end using a small coil producing a 
local magnetic field. This method works well in the numerical solution of the 
cold deflagration problem. However, such kind of experiment has not been 
performed yet. 
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Cold deflagration can be most likely observed on thinner crystals having 
a good thermal contact to the environment, so that the heat released in- 
side the crystal gets quickly removed and the temperature does not increase. 
As said above, the effect only exists within dipolar windows near tunneling 
resonances. 

It was shown that disorder in resonance fields of individual magnetic 
molecules is compensated for by adjustment of the dipolar field in the front, 
so that fronts of tunneling survive [30], 



1.3.4 Id theory of quantum deflagration 

Here we consider a more general situation in which the temperature of the 
crystal is increasing as the result of the decay of the metastable state, the 
case when the crystal is thermally insulated. The decay process is controlled 
by both the temperature (for any bias) and by the dipolar field (near tun- 
neling resonances) . The theory of the general quantum-thermal deflagration 
includes the relaxation equation ( |1.3[ ) and the heat conduction (energy diff- 



ision) equation (IL8|, as well as the expression for the dipolar field (1.51) in 



the Id approximation. The relaxation rate -T(T, B z ) was calculated for the 



generic Mni^Ac model ( 1.1 ) in Ref. [32] and for the realistic model of Mn^Ac 
containing the — AS Z term that splits tunneling resonances in Ref. |33|. 

Whereas an analytical solution of this problem has not been found, its 
qualitattive features can be well understood and the numerical solution based 
on discretization is available. In the case of a zero or weak transverse field, 
that was the case in all experiments up to date, spin tunneling is thermally 
assisted and it only modifies the main effect of regular deflagration, result- 
ing in tunneling peaks in the field dependence of the front speed v(B z ). As 
in the case of regular deflagration, ignition can be achieved by raising the 
temperature at an end of the crystal. 

Fig. 1 1 . 9 1 shows the front speed calculated for the bias and crystal size 
corresponding to the experiments in Refs. (IT -13 and using the relaxation 



rate r shown in Fig. 3 of Ref. 33 . The tunneling peaks arc quite pronounced, 
at variance with the results of these experiments. The latter can be due 
to a large ligand disorder in Mn^Ac that leads to a substantial scatter of 
the anisotropy constant D and thus of the positions of the resonances of 



individual molecules 61-63 , especially for the bias as strong as here. Just 
above 3 T and just below 3.5 T there are regions where the speed is too high 
to be measured in this calculation, an effect of ground-state tunneling. 

Spatial profiles of the magnetization, energy, and the total bias field in the 
deflagration front give an idea of the role played by spin tunneling. Fig. |1.10| 
shows the spatial profiles at the asymmetric peak of v at B z = 2.852 T in 
Fig. 1 1 . 9 1 Here the front speed is high because of tunneling at the face of the 



front, where in the lower panel the total bias field is flat at the level of the 
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Fig. 1.9 Numerically calculated speed of the deflagration front in a long Mni2Ac crystal 
for a weak transverse field. 





Fig. 1.10 Spatial profiles of the deflagration front in a small transverse field, B± = 0.05 
T at the peak of the front speed at B z = 2.852 T. There is a resonance spin tunneling at 
the face of the front and burning in its central and rear parts. In the tunneling region, the 
total field B Zj totsticks to its resonance value. 



tunneling resonance at £? Zj tot = 2.889 T. Magnetization distribution adjusts 
so that the dipolar field ensures resonance for a sizable group of spins that 
can tunnel. Tunneling of these spins results in energy release, the temperature 
and relaxation rate increase, and tunneling gives way to burning in the central 
and rear areas of the front. 

Formation of the asymmetric maxima of the front speed can be explained 
as follows. When B z increases, the peak of -Bz, tot that arizes due to the local 
dipolar field (central part of Fig. 1.6 1 reaches the resonant value. In thick 
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Fig. 1.11 Front speed for a strong transverse field (B± = 3.5 T) in the vicinity of the 
ground-state tunneling resonance at 0.522 T. The small peak on the left is due to the first- 
excited-state tunneling resonance. Left and right of the dipolar window the front speed is 
about 50 m/s 



crysta ls (R ^> Id) this happens if B z + B Z D = Bk, where B Z D is given by Eq. 
(1.531. This defines the left border of the dipolar window B z = Bk — B\ D 



(that differs from B z = Bk for the cold deflagration). At the left border 
of the dipolar window, a strong increase of v{B z ) begins. The maximum 
of i?z,tot sticks to the resonance value and becomes flat with progressively 
increasing width. Greater width of the resonance region results in a stronger 
tunneling and higher front speed. With further increase of B z , the right edge 
of the tunneling region moves too far away from the front core into the region 
where the temperature is too low. As the tunneling resonance in question is 
thermally assisted, it disappears at low temperatures, thus the flat region of 
-Bz,tot cannot spread too far to the right. As a result, the flat configuration 
of -E> 2 ,tot becomes unstable and suddenly B z-to t changes to the regular shape 
of the type shown in Fig. |1.6| 

If a strong transverse field is applied, the barrier becomes lower and it 
can completely disappear at a ground-state tunneling resonance. In this case 
r(T,Bk) is practically temperature independent and this maximum of the 
relaxation rate does not disappear at the highest temperatures achieved after 
burning, Tf. An example is the ground-state tunneling maximum at B z = 
0.522T in Fig. |1.5| Although at high temperatures this maximum is hardly 
visible in the log scale, it is clearly visible in the normal scale in Fig. 5 of 



Ref. 33 . In such strong transverse fields, the speed of the front becomes very 
high and spin tunneling plays the dominant role in the front propagation. 
Figure |1.11| shows a high front speed within a broad dipolar window 

B k - B^ <B Z < Bk + (1.62) 
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having the width of 125.5 mT. The front speed diverges towards the right edge 
of the dipolar window in accordance with Eq. (1.591 and becomes supersonic. 
A qualitatively similar behavior was observed earlier in calculations for the 
generic model of Mn^Ac, see Fig. 4 of Ref. 
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In contrast to thermally- 
assisted tunneling resonances, progressive flattening of B z ^ ot a,t its resonant 
value is not limited by the temperature before the front since ground-state 
tunneling occurs already at zero temperature. Thus the front speed diverges 
at the right edge of the dipolar window, Eq. (1.621, where the width of the 
tunneling region becomes very large. 

Comparing the present results with the analytical and numerical results 
for the cold deflagration, one can see that thermal burning in the central 
and rear parts of the front are stabilizing the process, so that the laminar 
solution, Eq. (1.59), holds up to the right edge of t he d ipolar window. There 

at Wo Xt 



1. 



1. 



is no breakdown of the laminar regime seen in Fig. 

Another feature of quantum deflagration is complete burning due to the 
temperature rise, in contrast to the incomplete burning in the cold deflagra- 
tion, Eq. (1.601. Although the speed of the cold deflagration front diverges at 
B z — > Bk + B Z D ^ (in the laminar regime), the amount of burned metastable 
population goes to zero, so that the rate of burning remains finite, Eq. (1.61 1. 
In quantum deflagration bur ning is complete [up to the equilibrium resudual 
population n( eq ) in Eq. Ql . 5p ] while the front speed is diverging, so that the 
rate of burning is diverging, too. 

Accordingly, the width of the front becomes very large at B z — > B^+B Z D \ 
in contrast to the width of the cold-deflagration front that remains constant. 
The structure of the front of the quantum-thermal deflagration near the right 
border of the dipolar window has a two-tier structure. First goes a fast front 
of tunneling that reverts a small fraction of the magnetization. The latter 
leads to heat release that ignites a front of thermal burning that burns all. In 
the stationary case the speed of the second part of the front is the same but 
it takes time to develop, thus the width of the whole two-tier front is large. 
Note that the speed of the quantum deflagration front is not limited by the 



speed of sound, contrary to the case of detonation 17 



1.3.5 3d theory of quantum deflagration 

As mentioned above, the Id theory of fronts of tunneling assumes a flat front 



that is not well justified because the dipolar field is given by Eq. ( 1.51 1 only 
at the symmetry axis. Different values of B Zito t away from the symmetry axis 
should self-consistently result in the distribution of the magnetization that 
depends on all coordinates x,y,z, i.e., in a non-flat front. 

On the top of this, there is an instability mechanism for a flat front at a 
smaller scale due to DDI. In Fig. |l.lO| we have seen that, approaching a front 
of tunneling from before, B z ^ot increases and reaches the resonance value, 
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Dipolar instability of a flat front of spin tunneling 




I Dipolar field 



Fig. 1.12 Dipolar instability of a flat front of spin tunneling. A leading part of the front 
(in the center) produces the dipolar fields on its neighbors that slow them down. 



then it becomes flat. Now. if a small fraction of the surface of a front (going 
from left to right and changing the magnetization in the positive direction) 
moves ahead of its neighbors, it produces a negative dipolar field on the lag- 
ging neigboring parts of the front, as any dipole, see Fig. |1.12| This brings 
the neighbors further from the resonance, so they tunnel later and their lag- 
ging increases. Conversely, lagging portions of the front produce a positive 
dipolar field on the leading part of the front that helps it to propagate faster. 
(The same mechanism leads to instability of flat domain walls considered in 
Ref. (29].) 

The DDI instability mechanism can potentially destroy any initially flat 
front of tunneling, making it microscopically rough. The question is whether 
micro-random dipolar fields produced by a micro-random magnetization in 
the front are still compatible with resonance tunneling. It is clear that rough- 
ness of the front breaks the concept of the adjustment of the system to the 
resonance, so that the speed of the front should decrease. On the other hand, 
spins are crossing the resonance, although at random times, so that still there 
should be a speed-up of the deflagration front near tunneling resonances. 

In 3d model of quantum deflagration the dipolar field was calculated using 



Eq. (1.541 for crystals of box shape with dimensions L x = L y -C L z using the 
relaxation rate r for B± = 3.5T shown in Fig. |1.5| The crystal was discretized 
with about 1 million total grid points in all 3 dimensions. The resulting system 
of first-order nonlinear equations was implemented in Wolfram Mathematica 
in a vectorized form using a compiled Butcher's 5th-order Runge-Kutta solver 
with a fixed step. 

As expected, roughness of the front due to the dipolar instability has been 



detected within the dipolar window, Eq. (1.62), where the computed front 
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Fig. 1.13 Front speed within the 3d model for a strong transverse field (B± = 3.5 T) in 
the vicinity of the ground-state tunneling resonance at B z =0.522 T. 




Fig. 1.14 Profile of the metastable population n in the 3d model of quantum deflagration 
for Mni 2 Ac at B± = 3.5T and B z = 0.5T (upper) and 0.56T (lower). 



speed is lower than within Id model, Fig. |1 . 11| Nevertheless, the front speedup 
due to spin tunneling is still huge, reaching sonic speeds in Mn^Ac on the 



right of the dipolar window, see Fig. 1.13 



Outside the dipolar window, a regular deflagration with a flat front and 
front speed v ~ 50m/s has been found for this value of the transverse field. 
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With entering the dipolar window from the left, the front becomes progres- 
sively non-flat with its central part leading. Front roughness emerges and in- 
creases with the bias. Fig. |1. 14] shows the profile of the metastable population 



n for the crystal with L x = L y = 0.2mm, as in experiments of Refs. 11-13 
for B z = 0.5T and 0.56T. The metastable population n is represented as a 
3d plot as a function of x and z with y — at some moment of time. The 
unburned cold portion of the crystal on the right is shown in blue, while the 
burned hot part on the left is shown in red. In the upper part of the figure 
showing the result for B z = 0.5T the front is essentially non-flat and there 
is some roughness, especially strong near the symmetry axis. The speed of 
this front v = 161m/s is already much greater than the speed of the regular 
deflagration, 50m/s. 

Numerical results for a larger bias B z = 0.56T and a longer crystal are 
shown in the lower part of Fig. |1.14| The front has a nearly sonic speed of 
v = 1674m/s and is very rough, while becoming flat again. The animation 
of this process looks like precipitation. Ignition of this front occurs at some 
distance from the left end of the crystal where the resonance condition is 
fulfilled. From this point, a very fast tunneling front is propagating to the 
right while a regular slow burning front is propagating to the left. 



1.4 Discussion 



Regular temperature-driven magnetic deflagration in long crystals of Mni2 
has been experimentally observed and is relatively well understood. The lack 
of a quantitative accordance between the theory and experiment can be at- 
tributed to still unknown temperature dependence of the thermal diffusivity 
k, as well as to the absence of a microscopic theory of relaxation in MM 
taking into account collective effects such as phon on/photon superradiance 
and phonon bottleneck. 

Effects of spin tunneling on ignition of deflagration and front speed near 
resonance values of the bias field have been experimentally detected in zero 
transverse field. However, these effect are due to thermally-assistent tunneling 
just below the top of the barrier and they are not strong. 

To the contrast, spin tunneling directly out of the metastable ground state 
in strong transverse fields can lead to huge effects such as supersonic quantum 
deflagration within the dipolar window around tunneling resonances. Unfor- 
tunately, creating an initial state for this process is practically difficult. In a 
strong transverse field also non-resonant spin tunneling is rather fast. During 
the system is being biased to reach the initial state close to the resonance, it 
is already relaxing and a large portion of the metastable population gets lost 
before a front of tunneling could start. In addition, non-resonant tunneling 
in a biased MM leads to heat release that can result in self-ignition if the 
crystal is thermally insulated. 
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It would be desirable to employ a fast field sweep to bring the MM into 
starting position for quantum deflagration in a strong transverse field without 
deteriorating its state. To observe non-thermal fronts of tunneling, thinner 
crystals with a good thermal contact to the environment have to be used. 
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